Glory Oscillations in the Index of Refraction for Matter- Waves 



(N 
O 

o 



in 



o 

O 

>V 

Oh: 



> 

OO 

o 
o 
m 
o 

(N 
O 

o 

>>: 
43 

Ph. 



- 1—1 

X 



Tony D. Roberts, Alexander D. Cronin, David A. Kokorowski, David E. Pritchard 
Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
(Dated: February 2, 2008) 

We have measured the index of refraction for sodium de Broglie waves in gases of Ar, Kr, Xe, 
and N2 over a wide range of sodium velocities. We observe glory oscillations — a velocity-dependent 
oscillation in the forward scattering amplitude. An atom interferometer was used to observe glory 
oscillations in the phase shift caused by the collision, which are larger than glory oscillations observed 
in the cross section. The glory oscillations depend sensitively on the shape of the interatomic 
potential, allowing us to discriminate among various predictions for these potentials, none of which 
completely agrees with our measurements. 
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named glory scattering by B. Cellini who observed the re- 
sulting circle of back-scattered light around the shadow 
of his head on the dewy grass. An even brighter glory can 
be observed surrounding the shadow of one's airplane on 
the clouds of spherical water droplets below Q. 

Glory scattering is not just an optical phenomenon. In- 
teratomic potentials are attractive at long range but have 
repulsive cores. Hence there is an impact parameter just 
beyond the potential minimum for which 0(b s i ory ) = 0. 
Quantum mechanics demands that atom-atom scatter- 
ing be treated as a wave phenomenon, but when the 
de Broglie wavelength is much smaller than the poten- 
tial, the semiclassical approximation can be used § in 
which the incident wave is viewed as traveling along clas- 
sical paths. There are then two possible paths to scatter 
into the forward direction — a long range diffractive com- 
ponent from large 6, and the aforementioned path near 
the potential minimum. Waves from these paths con- 
tribute to the scattering amplitude in the forward direc- 
tion, f(k,9 = 0), for which the partial wave treatment 
reads: 
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where the sums are both real. 

The behaviors of the real and imaginary parts of f(k, 0) 
are different. Im/ (the second sum in Eq. has a large 
value owing to the fact that sin 2 S averages to | where 
the scattering is strong (<$ ^> tt). The contributions from 
large / and l g i ory add on top of this large value, making 
Im/ less sensitive to the shape of the potential. On the 
other hand, the strong scattering region contributes zero 
on average to Re/, and the long range and glory contri- 
butions can either add or subtract, making it easier to 
discern effects from each region as well as their relative 
phase. 



The ratio, 
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has been shown to give information about the rate of 
increase of the interatomic potential V(r) for large r in- 
dependently of the strength of V(r) Also, oscilla- 
tions in p depend sensitively on the potential near the 
well bottom. This sensitivity has motivated several the- 
oretical predictions for p based on possible interatomic 
potentials |, |,|, §. 

We have measured p(k) in an atom interferometer us- 
ing a technique that parallels out previous measurement 
at fixed k |, |]. p is found by measuring Re/, which 
is related to the phase shift caused by the collision, as 
well as Im/, which is related via the optical theorem 
to the absorption cross section 0|. An atom interfer- 
ometer determines the phase shift and attenuation of 
the de Broglie wave in one arm of the interferometer 
that passes through a "target gas" by interfering with 
the other arm which does not. Specifically, propagation 
through a gas of length L modifies the wavefunction by 
a factor exp[i(n — l)ki a bL] where n is the complex index 
of refraction, n = 1 + 27rn<j/(A:, 0)/ki a bk, with rid the tar- 
get gas density and k and fc; a h the Na wavevectors in the 
center-of-mass and lab frames. 

The interferometer is composed of three nanofabri- 
catcd diffraction gratings forming a Mach-Zehnder ge- 
ometry |l(| . A gas cell exposes one path of the interfer- 
ometer to a gas while leaving the other path undisturbed 
(Fig. |). A beam velocity of 0.7-3.0 km/s (with -5% rms 
spread) is chosen by changing the carrier gas mixture in 
the supersonic oven. The velocity is measured to ±30 
m/s from the atom diffraction pattern of a grating. 

Several apparatus improvements were necessary in or- 
der to study the velocity-dependent index of refrac- 
tion. New 100 nm period gratings double the path 
separation — and hence the maximum usable velocity in 
the interferometer — compared to the old 200 nm grat- 
ings p2|. These gratings, developed by Tim Savas and 
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FIG. 1: Schematic of the atom interferometer and the gas 
cell. 



Hank Smith, have already lead to a number of advances 
in atomic and molecular physics |13j |. 

A thinner gas cell wall [|l4| also allows the atom paths 
to be closer together and hence allows higher velocities. 
The wall of the gas cell is a 10 fim thick Si wafer which 
is anodically bonded to a glass substrate. Channels cut 
into the glass create the volume of the cell and allow the 
beam to enter and exit (Fig. [l]) . 

The gas cell can be filled and emptied with a time 
constant of 1-2 sec (~10 times faster than before) us- 
ing computer-controlled valves and wider supply lines. 
Faster cycle times reduce errors due to phase drift caused 
by thermal and mechanical motion of the gratings. This 
drift has been further reduced by mounting the inter- 
ferometer on a rigid optical breadboard suspended by a 
vibration isolation system inside of the vacuum chamber. 

p is measured by comparing the interference pattern 
with and without gas in the cell. In the absence of gas, 
we observe an interference pattern in the detected signal 
Io that depends on the position x of one of the gratings, 



I {x) = N + A cos(^o + k g x), 



(3) 



where No is the average detected flux of atoms (ranging 
from 5-100 kHz depending on beam velocity), Aq is the 
amplitude of the interference pattern (typically 5-20% of 
No), and k g is the grating wavevector. <fo, the phase of 
the interference pattern in the absence of gas, is found by 
a fit to Iq(x) which also determines No and Aq. When the 
cell is filled with gas, the interference pattern becomes 

Igas(x) = N gas + A gas COs((f> gas + k g x) (4) 

from which p can be determined fet]: 
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A feature of this experimental method is that the 
measurements of attenuation and phase shift are rela- 
tive measurements, eliminating problems due to intensity 



fluctuations and phase drift. Furthermore, measuring the 
ratio of attenuation vs phase eliminates the dependence 
on the target gas pressure (which may fluctuate), and 
obviates the need to know it absolutely. It also cancels 
out the v~ 2 / b dependence of cross section and phase shift 
that can obscure the glory oscillations in /. 

In Fig. |[ p(v) is plotted as a function of Na beam ve- 
locity v for target gases of Ar, Kr, Xe, and N 2 . Each data 
point represents an average of 1-3 hours of data taken at 
one to three different pressures at a single beam velocity 
on a single night. Our previous index of refraction mea- 
surement j^] taken at a beam velocity of 1 km/s is shown 
for comparison. 

The improvements to the experimental apparatus men- 
tioned above have allowed us to acquire the large range 
of data in Fig. || with smaller statistical error than our 
previous measurements. We have also been more care- 
ful to consider possible systematic errors, correcting for 
these effects where feasible. 

The biggest source of systematic uncertainty is inter- 
ference from unwanted paths that reach the detector — 
for instance, the two paths that form the mirror image of 
the interferometer shown in Fig. [l], or from paths of Na 2 
molecules composing 25% of the intensity of the beam 
p2[ . To avoid molecule interference, we deliberately mis- 
align the gas cell wall at low velocities so that it appears 
wider than the Na 2 path separation. The 50 ^m-wide 
detector is positioned to avoid the other orders, but be- 
cause of beam alignment error and thermal drift, there 
is occasionally some detection of the wrong interference 
pattern, contributing as much as 5% to the interfering 
amplitude, Aq, except at higher velocities where it can 
be 50%. The resulting error in p is typically +3/— 0% 
but as much as +20/— 0% at the highest velocities of 
both the 100 and 200 nm data. This accounts for the 
asymmetrical systematic error bars in Fig ^[ 

Another systematic effect is the attenuation of the path 
outside the gas cell due to gas leaking from the ends of the 
cell. The integrated density of gas along the path inside 
the cell relative to the path outside is 120 ± 20 which 
results in a +1.7(0.3)% correction to the measured p. 

Another source of error arises from measuring the in- 
terference pattern while the gas cell pressure hasn't fully 
equilibrated after filling or evacuating. We do this to 
reduce the cycle time to ^20 sec leaving only 3 sec of 
dead time for filling or emptying. The resulting error in 
fitting to the time-averaged interference pattern, up to 
0.5% uncertainty in p, is included in the systematic error 
bars. 

Contaminants in the target gas must also be avoided. 
One such contaminant caused a 3% uncertainty in p for 
the 200 nm Ar data before it was fixed. 

Error bars on the new data points in Fig. ^| include 
the sum of statistical and systematic errors. Because the 
only non-negligible systematic comes from beam align- 
ment and drift, we have an independent check of this 
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FIG. 2: p as measured for Na waves in Ar, Kr, Xe, and 
N2 (• using 200 nm gratings, o 100 nm), showing evidence of 
glory oscillations in comparison to p as derived from predicted 

potentials: Na-Ar Q{— ), &L ), §(■ • ■ ), {jjL-r ")' 

[ LSI C ); Na-Kfll(— )7||( ); and Na-Xe @](— ), 

[19|( ), pl|(- ■ ■ ). Our 1995 measurement is also shown 

[ |f(boxes) . The horizontal tick marks show the statistical er- 
ror. The vertical error bar line includes systematic errors. 



error by examining the repeatability of data taken on 
different nights. We expect that the error in alignment 
or drift may exceed our estimates on some nights, but 
this can be recognized: p will be artificially low com- 
pared to other nights, we will see this deviation in all the 
gases measured, and the measurement will have larger 
statistical error bars as the ratio changes due to drift. 
An example of these features can be seen in the severely 
outlying data points measured at 1450 m/s for Kr, Xe, 
and N2. 

Another check of systematics can be made by compar- 
ing the 100 and 200 nm data where they overlap. At this 
velocity, unwanted interference is seven times larger for 
the 200 nm data than for the 100 nm data, yet the agree- 
ment between the two is well within the error estimates. 

The data are also consistent with older measurements 
taken at v = 1 km/s ||, with the exception of Xe, which 
disagrees dramatically. However, no care was taken in 
that experiment to consider systematic errors, which we 
estimate were +30/— 10% in p due to the unwanted in- 
terference of Na2 molecules and could have been worse 
for this particular measurement if thermal drift was ex- 
cessive. 

To compare our measurements with a prediction for 
V(r), the sum in Eq. ([j]) must be computed. The Na 
wavelength is much smaller than the range of the in- 
teratomic potentials we are considering, so hundreds of 
partial waves contribute to the sum and we are justified 
in replacing it with an integral: 



Re/(fc) = k b sin 28(b) db 
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Im/(fc) = k 26 sin 2 8(b) db 
Jo 



(6a) 



(6b) 



where b = (I + ^)/k is the classical impact parameter. 

Because the kinetic energy is much greater than the 
well depth of the interatomic potential, we can find 8(b) 
using the Eikonal approximation — that 8 is the accumu- 
lated phase shift of the atom traveling at constant speed 
along a straight-line path with impact parameter b: 



8(b) 



2kh 2 



dzV(Vb 2 + z 2 ), 



(7) 



where p, is the reduced mass. Predictions using this ap- 
proach are valid to 6% in p for the Na-Ar system and 3% 
for the Na-Kr and Na-Xe systems in comparison to an ex- 
act quantum treatment p3| . Predictions for p must also 
average over the room-temperature distribution of target 
gas velocities in our experiment, which damps the glory 
oscillations somewhat at lower beam velocities lj, [?]] . 

Fig. U also shows calculations of p(v) based on pre- 
dictions of V(r) for Na-Ar R ^ @> Na " Kr 
|5|, [19), and Na-Xe |§ fH], ^[derived from spectro- 
scopic measurements and beam scattering experiments. 
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The oscillations in these curves are glory oscillations, and 
are evident in the data. Even after allowing for statisti- 
cal and systematic errors, our data disagree with these 
predictions. We believe this is because our measurements 
are very sensitive to the shape of the potential near the 
minimum (see examples in ^4|), where the transition 
from the repulsive core to the Van der Waals potential is 
poorly understood, even though the well depth and loca- 
tion of the minimum are fairly well known. The N2 data 
are remarkable as well for showing a large variation in p 
without any evidence of an oscillation, unfortunately, the 
potential is not understood in enough detail to explain 
this behavior |26| . 

Because there is good agreement on the location and 
depth of the potential minima, we can use this knowledge 
to scale our measurements of p in a way that is sensitive 
to the shapes of all three Na-rare-gas potentials, inde- 
pendent of the well location and depth. We can write 
any potential as V(r) = D e g(r/r e ) where D e and r e are 
the depth and location of the minimum and g(x) is a 
dimensionless function with a minimum g(l) = — 1. In 
the Eikonal approximation p depends only on g(x) and 
the dimensionless parameter a = D e r e /hv. a is approxi- 
mately the phase shift in radians accumulated along the 
glory path. Changing the well depth or location of a po- 
tential without changing its shape g(x) is equivalent to 
changing the velocity scaling of p(v). 

In Fig. H we have plotted measurements of p for Ar, 
Kr, and Xe verses the dimensionless parameter a, using 
values of D e and r e taken from Because of thermal 
averaging, the values of p we measure also depend weakly 
on the ratio of v to the mean speed of the scatterer, but 
differences in p(a) due to thermal averaging would be at 
most of 5%. Also included are data for Ne, with new sys- 
tematic error estimates, which were taken with the older 
apparatus j§ §. Plotted in this manner, the combined 
data clearly show a full glory oscillation, as well as a hint 
of a second. For comparison we show p as derived from 
the Lennard- Jones potential, g(x) — x~ 12 — 2x~ 6 , and 
the Morse potential g[x) = e -2 ^ -1 -' — 2e~^ x ~^ using 
(3 = 4.7 which is appropriate for the Na-Ar system |Q. 
The continuity of the data plotted together in this fash- 
ion demonstrates the similarity of Na-rare-gas potentials 
for these four systems and shows that Lennard-Jones and 
Morse potentials are too simplistic to accurately model 
the data here. 

In conclusion, we have measured p to a typical accu- 
racy of 9% at various Na velocities in gases of Ar, Kr, Xe, 
and Nj. The consistency of our scaled data for different 
gases suggests that the shape g(x) of different Na-rare- 
gas potentials is similar. However, discrepancies in the 
comparison to predictions of p suggest that the shape of 
the potential is not well understood theoretically. Unfor- 
tunately, the difficulty of the inverse scattering problem 
p7| , p8| prevents us from deriving the potential based on 
knowledge of f(k). However, our explorations of modifi- 
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FIG. 3: p for the Na-rare-gas systems [Na-Ne (boxes), Na- 
Ar (bullets), Na-Kr (diamonds), and Na-Xe (black triangles)] 
plotted versus the dimensionless parameter a = D e r e /hv. 
Also shown is the predicted p assuming a Lennard-Jones po- 
tential ( — ), and a Morse potential ( ). 



cations to suggested interatomic potentials do indicate 
that they would agree better with our measurements 
if V(r) were modified such that V(ro) = occurs at 
smaller tq and the well were made less deep in the range 
r e < r < 2r e [pij . Our measurements should be useful in 
evaluating future refinements of the relevant interatomic 
potentials. 
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